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Abstract 

To a first approximation, the quark mixing matrix has Of^ = = Oi whereas the 
lepton mixing matrix has ^23 = tt/A. We show how this structure may be understood if 
the family symmetry is Qs, the quaternion group of eight elements. We find three viable 
scenarios for the Majorana neutrino mass matrix, each depending on 4 parameters and 
predicting a specific mass spectrum. The phenomenology of the two Higgs doublets 
which generate the Yukawa sector is analyzed and testable predictions are derived. We 
discuss also the closely related model based on L'4, the symmetry group of the square. 
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There are 3 families of quarks and leptons with accompanying 3x3 mixing matrices: 
VcKM = Vj^Vd hnking the {u,c,t) quarks to the {d,s,b) quarks, and Umnsp = UiUy hnk- 
ing the (e, /i,r) charged leptons to the (z^i, z^2, t's) neutrinos. The Vckm matrix may be 
parametrized with 3 angles 9^^ and 1 phase 5' and similarly Umnsp, after absorbing two 
relative Majorana phases in the neutrino mass eigenvalues. Numerically, \Vus\ — 0.22, 
\Vcb\ — 0.04, \Vub\ — 0.004, thus setting 9% = 6'23 = is a good first approximation of 
Vckm- On the other hand, Umnsp has a very different structure: both 9[2 and 6*23 are 
known to be large, with = as its best value experimentally (and 9{2 7^ 7r/4). 

We propose in this letter a new understanding of Vckm and Umnsp in terms of the non- 
Abelian discrete symmetry Qs, the quaternion group of 8 elements. Using a natural generic 
assignment of quarks and leptons and the simplest nontrivial Higgs content, we show that 
there are only 4 possible scenarios for leptons (but only 3 are viable phenomenologically) and 
1 scenario for quarks. We will also discuss the closely related discrete symmetry D4, which 
is the symmetry group of the square. The latter also has 8 elements and the same character 
table as Qs with 5 conjugacy classes and 5 irreducible representations (irreps): 

1 , and 2 (see Table Q). Two specific models based on D4 x Z2 have recently been 
proposed [H |2] . 

Table 1: Character table of Qs (-D4). Here n is the number of elements in each conjugacy 
class, while h is the order of any element g in that class, i.e. the smallest integer such that 
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Before we show how it all works, let us present our main results. Using Qs, we find 
the following 4 scenarios for the Majorana neutrino mass matrix (in the basis where the 
charged- lepton mass matrix is diagonal): 
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Scenarios (1) and (4) correspond to mass matrices with two texture zeros, which have been 
studied in the hterature 0; in particular scenario (4) is known to be ruled out. Scenarios 
(2) and (3) both have = "^33 and one texture zero. These are new structures which are 
also viable phenomenologically as we will show. 

The group Qs may be generated by the eight 2x2 matrices ±1, iiai, ±icr2, ^icr-s, 
whereas D4 may be obtained with ±icri^3 replaced by ±(Ti^3. Each set is also a faithful two- 
dimensional irrep 2 of the respective group. Geometrically, the group Qs may be associated 
with the 8 vertices of the hyperoctahedron (dual of the hypercube) in four dimensions. 
The two superscript signs for the one-dimensional irreps correspond to the characters of 
the C4 and C5 classes, generated by iiai ^ (or ±o"i^3) respectively. Among themselves, 
the 4 one-dimensional representations transform as Z2 x Z2. In the case of Qs, the irreps 
l'^ , 1"^, and 1 are completely interchangeable, since they share exactly the same group 
properties. This means that if a theory contains a set of these irreps, replacing them with 
those obtained by any {S3) permutation of the 3 classes 6*3,4,5 will not change the physical 
predictions of the theory. Specific examples will be given below. In the case of D4, only 
and l""*" are equivalent, since the conjugacy class C3 is distinguished from C4 and C5 by 
the order of their elements (see Table H}. For both groups, the basic tensor product rule is 
2x2 = l"*""*" -|- l"*"" + l""*" + 1 , but the doublet components are combined in different ways 
for Qs and D4, as shown in Table |21 Since the two-dimensional irrep of D4 is real, ("i/'i, ^'2) 
transforming as a doublet implies that {ipl, ^/'D is also a doublet, whereas in the case of Qs, 
the correct assignment is {ip2y expected, just like a doublet under SU{2). 

Under Qs, we assign the 3 quark and lepton families and two Higgs doublets as follows: 
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Table 2: 2 x 2 decompositions of Qs and D^. 
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As a result, each quark mass matrix in the basis (1 , 1" 

fa d 0\ 



is of the form 
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where a, 6, c are proportional to (05) and e to (02)- This means that the third family does 



not mix with the other two, i.e. 



0, which is a good first approximation. On the 
other hand, the charged-lepton mass matrix in the basis (1''"'', 2) is given by 

/a \ 

X; = c 6 , (9) 

vO -h -CI 

where a, 6 are proportional to (05) and c to (02)- This matrix is easily diagonalized by a 
rotation of 7r/4 on the left and on the right, i.e. /x, t = (/2 ±/3)/v^ and if ^ = (^2 + ^3) / 



with rrip 



b\, and nir = \c + b\. 



The neutrino mass matrix is assumed to be Majorana and generated by the naturally 
small vacuum expectation values (VEVs) of heavy Higgs triplets = , ^ ^i) ■ The 
assignment of triplets to the different Qs irreps is crucial for the resulting neutrino mass 
pattern. 

Scenario (1). Since the triplet VEVs are induced ^ via trilinear couplings of the 
form ^i0j0fc, the requirement of Qs symmetry would allow only ^1 ~ 1++ and ^2 ~ 1^ 
to contribute. In that case, Aijy has \m2\ = jmsl and z/g is unmixed, which is not realistic. 
However, Qs is expected to be broken softly, thus the scalar trilinear C,i(j)j4'k terms may induce 
small VEVs also on (1^3, .^4) ~ 2, by which 
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where a comes from (^°), b from (^2)? ^ = ^(^4); and / = —h{C,^). In the e,fi,T basis with 
l-i'^'T = (h^h) I V^; the neutrino mass matrix takes the form of Eq. ((H), with c,d = (e±/) / \/2. 
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Two texture zeros in the t'^z/^ and Vt-Vt entries are thus derived (for the first time) by our 
apphcation of the Qg family symmetry. This is known to be a viable pattern jSl E| and in 
most cases predicts an inverted mass spectrum. Here = tt/A and = are obtained 
in the limit of c = ±d. Deviations from maximal 2 — 3 mixing are allowed proportionally 
to the nonzero value of ^^3, which can be as large as the experimental upper bound. With 
the present experimental constraints, we find \m2\ > 0.04 eV and \m^\ > 0.015 eV, as shown 
in Fig. 1. The neutrinoless 2/?-decay rate is controlled by niee = \a\ > 0.02 eV. A quasi- 
degenerate spectrum can be obtained, when in Eq.(^ a ^ b and c,d are much smaller j^; 
in this limit the ordering of the spectrum can be also normal (see Fig. P). 
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Figure 1: The allowed region in m2 — plane for scenario (1) is presented in eV units. 
The masses are scanned in the experimental allowed range : Am^^i = (7.7 — 8.8) x 10^^ eV^, 
Aml^^ = (1.5 - 3.4) x lO^^ eV^ tan^ ^^^i = 0.33 - 0.49, sin^ 2^atm > 0.92, and sin^cHooz < 
0.2. 

There are 3 equivalent assignments of Higgs doublets and triplets which result in scenario 
(1), as listed in Table El In the choice B2, J^i is diagonal if we redefine /2 3 as 1^ 2^ 
has the form of Eq. (pj) automatically as it should. In the choice C3, the charged-lepton and 



5 



Table 3: Assignments of the Higgs doublet 02 and of the Higgs triplets and ^2 corre- 
sponding to the four scenarios for the neutrino mass matrix described by Eqs. (0)-©- The 
assignments (pi ~ and {(,3, (,4) ~ 2 are kept fixed in all cases under study. 
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neutrino mass matrices are given by 
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Diagonalizing Aii by 

/I \ 

U= -i/V2 I/V2 , (12) 
VO t/V2 I/V2) 

we find that A^(f'^'^) = UMuW^ has the same form as Eq. dTj), with c,d = {f ± ie)/\^. In 
the quark sector, if we replace the 3 one-dimensional irreps of Eq. ((H)) with any other 3, and 
put the remaining one in the lepton sector, we again have the same physical predictions, i.e. 
Eqs. © and (P). 

Scenario (2). Instead of using Higgs triplets of the same one-dimensional irreps as the 
Higgs doublets which we did in scenario (1), consider the use of one triplet in the 1 irrep 
(as 01 ) and a second triplet in one irrep equivalent but different from that of 02. There are 6 
such assignments, as shown in Table 01 which are all equivalent as expected. For definiteness, 
let us study the choice Bl. In this case, Aii is given by Eq. © and Aiu by Eq. ((T)), so that 
in the e, /i, r basis we obtain Eq. (0) (with a redefinition of c and d). Let us now prove that 
this new pattern is viable phenomenologically. If c and d are relatively real, than m^, = b and 
^3 = S23i^^ + C23VT, with S23 = —d/ycF+d?, and C23 = c/ ycF+d? , i.e. O^^, = necessarily 
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and tan = —d/c is arbitrary (the second x Z2 model j2] also has this property, but its 
other predictions are different). In this limiting case we have the identity a = mi + m2 — 
and the sum rule sfg^^i + £^2^2 = ""^3, which implies inverted hierarchy as shown in Fig. 2. 
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Figure 2: The allowed region in 1712 — plane for scenario (2) is presented in eV units, 
in the limit = 0. The masses are scanned in the same experimental allowed range of 
neutrino parameters as Fig. ^ 



If c and d are not relatively real, ^^^3 becomes nonzero, and the resulting allowed region 
in the m2 — plane is similar to that of Fig. ^ In particular, a degenerate spectrum is 
allowed, when the matrix (jSJ becomes close to the identity matrix (a ~ 6 and c, d much 
smaller) [S]. 

Scenario (3). The last viable scenario under Qs has one Higgs triplet in the same irrep 
as 02 and a second triplet in one of the other two equivalent irreps. There are 6 equivalent 
assignments, as listed in Table El Consider the choice El. In this case, A4i is given by Eq. Q 
as before, but the neutrino mass matrix is now 
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Since this is the same as Eq. ()1U|) except for the placement of a, it becomes the matrix of 
Eq. Q in the e,/i, r basis, in analogy with Eq. (^. This is another new viable pattern and 
predicts a normal hierarchy. Here ^23 = and = are again obtained in the limit 
c = ±d. With the form of Eq. (jHl), the constraint sin6'i3 < 0.2 implies that sin^ 26*23 > 0.987, 
and predicts 



Consider the choice F3. In this case, Aii is the one given in Eq. (fTTjl . which is diagonalized 
by U of Eq. (fT^ . On the other hand, the neutrino mass matrix is of the form of Eq. (jH)) 
already, so isn't the extra rotation of 7r/4 required by U going to change its form? The 
answer is no, because 



which also rotates the — 1/^ sector by 7r/4. As a result, Eq. Q retains its form: the 
observed maximal atmospheric mixing emerges from maximal mixing in both neutrino and 
charged lepton sectors! It is easy to work out the details in each of the other cases as well. 

In all 3 scenarios, i.e. Eqs. (P), and Q, there are 4 parameters in Ai^'^'^\ This 
means that the 3 neutrino masses and their 3 mixing angles are related. In particular, the 
absolute scale of mi.2 3 is constrained in each shown in Figs. 1,2, and Eq. (jl4|) . Thus 

they are all testable predictions. 

So far we have a successful Umnsp matrix but only an approximate Vckm matrix because 
^13 = ^23 = 0- This means that the quark sector has an extra global symmetry associated 
with the third family, i.e. a top/bottom flavor number. This symmetry may be broken, for 
example, by adding to the specific assignment of Eq. ((Tj) the Higgs doublet 03 ~ 1 which 
contributes to the (23) and (32) entries of A^g in Eq. (jS)) to allow ^23 and also Of^ 7^ 0. 
This will of course also affect Aii of Eq. Q but only in the fi — t sector. It means that 
even if 6*^3 = in scenarios (1) and (3), ^33 should deviate from 7r/4 as well, having the same 
origin as deviations of ^23 and 6*^3 from zero. (For a possible indication of deviation from 
maximal atmospheric mixing in SuperKamiokande data see for example Ref. [0].) 

If we use D4 instead of Qs, then according to the multiplication rules of Table 2, we find 
one realization of scenario (1), i.e. the choice C3 of Table El two of scenario (2), i.e. the 
choices CI and C2, and two of scenario (3), i.e. the choices Dl and D2. There are also two 
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other viable scenarios (but with 5 parameters), i.e. 
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Scenario (3') is obtained with the choices A2,A3,B1,B3 of Table El and scenario (4') with 
A1,B2. If we assume c = d in scenario (3'), then we obtain the result of the first x Z2 
model P, but without using the extra Z2. 

In the Higgs sector, 0i and 02 are distinguished by an odd-even parity. The ensuing 
scalar potential is well-known, having in general a minimum with nonzero VEVs for both 05 
and 02. Consider A^g of Eq. (jH}. The 2x2 sub-matrix spanning the first two families may 
be diagonalized in general by a rotation on the left and a rotation on the right. Specifically, 
in the case of the d and s quarks, it may be written as 
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Since 0? couples to the diagonal entries and 0° couples to the off-diagonal entries, there are 
ds'^ and sd'^ couplings to each given by 



. D {[-SrCr{cI- sl)md + SlCl{cr- s%)ms] ds"" 
+ [-slCl{c% - s\)md + srCr{c\ - s\)ms\ sd"^ 
X [<PVv2 - 0?/^l) + H.c. 

Note that, even though M.i of Eq. Q is not diagonal, the /x — r sector has c\ = s| 
s|. = 1/2. Hence Eq. p7|) shows that flavor-changing fi — t interactions are absent. 



(17) 



In our model, the state (fi0? + V24>2) / y vl + V2 is identifiable with the neutral Higgs 



boson of the Standard Model. Its orthogonal state = (fi02 — V2(pi) / \J vf + V2 appears 
in Eq. (fT7|) and couples to both ds^ and sd^, thereby contributing to the Ki — Ks mass 
difference. This contribution depends on the unknown parameters sr and ruh. It is suitably 
suppressed if either sr is small or m/j is large. For example, if sr is negligible, then [7] 
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123 GeV, si ~ md/ms, Bk = 0.4, Jk = 114 MeV and = 7 MeV, 



this contribution is 1.1 x 10 ^^(100 GeV/m/i)^, the experimental value being 7.0 x 10 
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Similarly, the contribution to /^nici/mj:, is estimated to be 1.3 x 10~^^(100 GeV/m/i)^, which 
is well below the experimental upper bound of 2.5 x 10~^^. In the leptonic sector, is 
predicted to have the interaction 

+ + t) ("^M^^'^ + ^r/X/x'^)} + H.C. 

Its decay into t^t^ and pairs will be crucial in the verification of this model. 
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